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ABSTRACT 

We report the discovery of a protocluster at z ~ 6 containing at least eight cluster member galax- 
ies with spectroscopic confirmations in the wide-field image of the Subaru Deep Field (SDF). The 
overdensity of the protocluster is significant at the 6cr level, based on the surface number density of 
i'-dropout galaxies. The overdense region covers ~ 6' x 6'(14Mpc x 14Mpc in comoving units at 
z = 6), and includes 30 i'-dropout galaxies. Follow-up spectroscopy revealed that 15 of these are real 
z ~ 6 galaxies (5.7 < z < 6.3). Eight of the 15 are clustering in a narrow redshift range (Az < 0.05 
centered at z = 6.01), corresponding to a seven-fold increase in number density over the average in 
redshift space. We found no significant difference in the observed properties, such as Lya luminosities 
and UV continuum magnitudes, between the eight protocluster members and the seven non-members. 
The velocity dispersion of the eight protocluster members is 647 ± 124 km s^ 1 , which is about three 
times higher than that predicted by the standard cold dark matter model. This discrepancy could 
be attributed to the distinguishing three-dimensional distribution of the eight protocluster members. 
We discuss two possible explanations for this discrepancy: either the protocluster is already mature, 
with old galaxies at the center, or it is still immature and composed of three subgroups merging to 
become a larger cluster. In either case, this concentration of z = 6.01 galaxies in the SDF may be one 
of the first sites of formation of a galaxy cluster in the universe. 

Subject headings: early universe — large-scale structure of universe — galaxies: high-redshift - 
galaxies: clusters: general 



1. INTRODUCTION 

Exploring the structure formation and evolutionary 
history of the early universe is an issue of strong cur- 
rent interest in astronomy. In the cold dark matter 
(CDM) model, clusters of galaxies form in the densest 
peaks of dark matter in the early universe. These grow 
by merging and by accr eting material from surround- 
ing low-density regions (Sprin gel et al.l [2005). There- 
fore, the growth of density perturbations measured by 
the comparisons of abundances and mass distributions of 
present-day clusters with those at earlier times provides 
uniqu e constraints on the ACDM concordance model 
fe-g UVikhlinin et al.l [20031 IVoul [20051 : iMortonson eTaLI 
120 111) . Theoretical models predict that galaxies lying 
inside these high-density regions may have formed ear- 
lier and/or evolved more rapidly, compared with their 
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surroundings (e.g., i Kauffmann 1995; Benson et al.H200lT ; 
iThomas et all 120051: iDe Lucia et al.l 12006( 1 . Clusters of 
galaxies are, thus, the noteworthy site of dark matter 
structure formation and galaxy evolution. In the lo- 
cal universe, cluster galaxies represent a distinct rela- 
tionship in the color-magnitude diagram. The "red se- 
quence" in clusters is composed mainly of spheroidal 
and lenticular galaxies w ith old stellar population s and 
high stellar masses (e . g.. IVi svanathan & Sandagd ll977t 
iGladders fc Yei I2000J: iLerchster et all 120111) . Further- 
more, brighter and massive galaxies are likely to lie the 
central region of a cluster. When and how, in the history 
of the universe, were these distinct properties formed? 
Highly evolved clusters have been intricately affected by 
both nature and nurture; thus, it is difficult to sort out 
which factors are most important to cluster and galaxy 
evolution. One promising line of study is to investigate 
directly the primitive properties of galaxies belonging to 
the first clusters of galaxies in the early universe. Proto- 
clusters in the early universe would provide a great deal 
of information on the primordial conditions of clusters at 
their birth. 

Star-forming galaxies, such as Lya emitters (LAEs) 
and Lyman break galaxies (LBGs), are almost the only 
tracers that have been used to follow the evolution of 
large-scale structures at high redshifts. We have very 
few clear examples of these star- forming galaxies strong! 
clustered beyond z — 3 (e.g., iLe Fevre et al.l 1199 
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2009( 1. Subsequent se arches have been carried out at even - 
higher redshift s fe.g , IVenemans et alJl2002l 12 004. 2007: 
IQverzier et all 120061 120081 ICapak et all 120111) Over- 
dense regions of i'-dropout galaxies out to z ~ 6 have 
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been found tiMalhotra et all 120051: iStiavelli et~al] 120051: 
iZheng et al] 120061: iKim et a l. 2009); however, most of 
these potential protocluster members have not been iden- 
tified through spectroscopy. Analysis based only on 
photometric data limits a detailed understanding of the 
physical processes, such as three-dimensional structure, 
in the high-density environment. Furthermore, sample 
could be contaminated by lower-redshift objects; it di- 
lutes intrinsic properties of protoclusters. 

When cluster formation approaches its final stage, 
galaxy evolution is predominantly influenced by envi- 
ronmental effects. Some differences between the proper- 
ties of protocluster galaxies and field galaxies appear at 
z ~ 2—3: the masses o f cluster galaxies are higher than 
those of field galaxies dSteidel et ahl 120051 : iKuiper et al.l 
2010; Hatch et al. 2011), and the fraction of galaxies har- 
boring an active galactic nucleus (AGN) i n protoclusters 
at z ~ 3 is higher than that in the field (jLehmer et al.1 
2009). These findings indicate that the first step in 
forming a cluster takes place at a much higher r edshift . 
In earlier epochs, at z ~ 4—5, lOverzier et al.l (|2009b) 
found no significant difference in the stellar mass be- 
twee n protocluster g alaxies and field galaxies. Addition- 
ally, lOverzier et al.l (|2009bT ) showed that the total stellar 
mass of protocluster galaxies at z ~ 4 is much smaller 
than that of massive cluster galaxies at z ~ 1 which are 
expected to have formed at z > 4. This may suggest that 
there are more massive structures at z > 4, which will 
become rich cluster in the local universe. It is impor- 
tant to investigate protoclusters at even higher redshifts 
directly, to understand how the clusters are formed and 
how galaxy evolution depends on the environment. In- 
vestigating protoclusters during the epoch of reionization 
may also be important for understanding the early his- 
tory of the universe, because the reionization process and 
thermal histories in overdense regions are expected to be 
remarkably different from those of average density, due 
to the enhanced numbe r of ionizing sources and higher 
radia tive feedback (e.g. jMcQuinn et al. 2007; Iliev et aLI 
12001) . 

Most of these protoclusters have been discovered in re- 
;ions centered on radio gal axies (RGs) or quasars (QSOs) 
Milev fc De Breuckll2008l ). RGs and QSOs have been 
used as potentially useful probes of large-scale struc- 
ture, as signposts of possible regions of galaxy over- 
density. However, strong radiation from RGs or QSOs 
may provide contradictory feedback that suppresses 
nearby galaxy formatio n, especially of low-mass g alaxies 
(|Barkana fc Loeblll999D . iKashikawa et al] (|2007li found 
a ring-like structure of LAEs around a luminous QSO at 
z ~ 5, possible evidence of negative feedback on star- 
formation activity in galaxies neighboring QSOs. Due 
to their low number density, it is difficult to discover 
protoclusters without RGs or QSOs at z > 3. A few pro- 
toclusters or large-scale structures have been discovered 
serendipitously in rand om fields at high redshift (e.g., 
Steidel et ail Il9 98: Shi masaku et al.l 120031: lOuchi et al] 
20051 : iLemaux et aljl2009t ). suggesting that early massive 



structures do not always host RGs/QSOs. 

In this paper, we present the most distant protoclus- 
ter candidate known, at z = 6, hosting at least eight 
spectroscopically confirmed members. This protocluster 
was found from an examination of data taken with a 
wide field-of-view (FOV) camera mounted at the prime 



focus of the Subaru telescope and aimed at the Subaru 
Deep Fi eld (SDF; RA=1 3 h 24 m 38!9, DEC=+27°29'25'.'9 
[J2000]: IKashikawa et all (120041) V The pri me focus cam- 
era, "SuprimeCam" (jMivazaki et al.l 120021 ) . has an FOV 
of 34' x 27', and is a powerful instrument for discover- 
ing similar rare objects. This paper is organized as fol- 
lows. §2 describes the imaging data used in this study 
and our selection of z ~ 6 galaxies. In §3, we deter- 
mine the significance of the surface number density of 
the z ~ 6 galaxies in the overdense region. In §4, we 
describe our follow-up spectroscopic observations. In §5, 
we discuss the properties and the structure of the pro- 
tocluster. The conclusions are given in §6. We assume 
the following cosmological parameters: fijyr = 0.3, Qa = 
0.7, Ho = 70 km s _1 Mpc -1 , which yield an age of the uni- 
verse of 910 Myr and a spatial scale of 40 kpc arcsec -1 in 
comoving units at z = 6. Unless otherwise noted, we 
used comoving units throughout. Magnitudes are given 
in the AB system, and we used a 2" aperture. 

2. SAMPLE SELECTION 
2.1. Photometric Data 

We used the SDF public data, with limiting magni- 
tudes of B = 28.45, V = 27.74, R = 27.80, i' = 27 A3, 
and z' = 26.62 (2" aperture, 3a). We have also obtained 
new, deeper R-, i'-, and z'-band images. These images 
were constructed by stacking all the data taken from 2001 
to 2008 in the course of a study of dist ant supernovae 
(jPoznanski et al.1120071 : iGraur et al.ll2011l ). containing al- 
most 30 hours worth of integration time in total. The 3a 
limiting magnitudes of these new deep images are 28.35, 
27.72, and 27.09 at R-, i'-, and z'-bands, respectively. 
These are about 0.5 mag deeper than the SDF public 
data. We selected z'-dropout objects using these deep 
R-, i'-, z'-band images and the public B- and U-band 
images. All five images were convolved to a common 
seeing size of C'98. We obtained J-b and imaging of the 
SDF with the WFCAM on UKIRT (jCasali et al] 12001 
in March, April, and July 2010 (Hayashi et al., in prep). 
The SDF was entirely covered by a mosaic J-band image, 
although the depth was not uniform from field to field. 
We supplemented our data with this J-band imaging, the 
limiting magnitude of which was > 2 mag shallower than 
the z'-band image. The seeing size is l'.'l. This allowed 
us to reject apparent contamination by M/L/T dwarfs. 
The pixel scales of all six images are 0'.'20, and the details 
of the image properties are summarized in Table [TJ 

We performed object detect ion and photometry by 
running SExtractor (ver. 2.5.0) ([Bertin fc Arnoutill996l) 
on the images. Object detections were made in the z'- 
band. Then, the magnitudes, and several other photo- 
metric parameters, were measured in the other bands at 
exactly the same positions and with the same apertures 
as in the detection-band image. This task used the SEx- 
tractor "double image mode." We marked objects that 
had five connected pixels with signal that was 2a above 
the sky background RMS noise. Photometric measure- 
ments were made at the 2a level. Objects detected in 
regions with low S/N ratio, near the frame edges or near 
saturated pixels around bright stars, were removed from 
the catalog. The remaining effective area of analysis was 
876 arcmin 2 . Finally, ~ 102000 objects were detected 
down to z' = 27.09 (3a limiting magnitude). To esti- 
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TABLE 1 
Photometric Data 










B 


V 


R 


i' 


z' 


J a 


^Zim,3cr 

integration time (hour) 


28.57 
9.9 


27.85 
5.7 


28.35 
27 


27.72 

27.7 


27.09 
30.9 


23.30-24.80 
1.1-10 



a J-band image is not uniform over the field, due to mosaicing four regions. The protocluster region is located in the deepest portion of the J-band 
image. 

3ct limiting magnitude in a 2" aperture 



mate the detection completeness of the z'-band image, 
we used the IRAF task mkobjects to create artificial ob- 
jects on the original image. Artificial objects were cre- 
ated with the same FWHM (0'.'98) as real images, and 
were randomly distributed. To avoid blending artificial 
objects with real objects, we avoided positions close to 
the real objects with distances shorter than 1.5 times the 
FWHM of the real objects. We extracted the artificial 
objects using SExtractor with the same parameter set. 
We generated 3000 artificial objects in the 23 to 29 mag- 
nitude range, and repeated this procedure 20 times. The 
detection completeness was more than 90% at z' = 25 
and 70% at z' — 27. 1 (3cr limiting magnitude). 

2.2. Selection Criteria of z ~ 6 Galaxies 

Figure Q] presents a color-magnitude diagram of the 
z'-band detected objects. We selected z ~ 6 galaxy can- 
didates using the following criteria: 1) i' — z' > 1.50 
down to z' < 27.09 (The la limiting magnitude of i'- 
band was used as i' — z' color limit, if objects were not 
detected in i'-band image), 2) null detection (< 2a) in 
B- and V^-bands, 3) z' — J < 0.80 (this criterion was 
imposed only for objects detected in J-band.), and 4) 
z' > 25.0. The 1 st color criterion picked up galaxies with 
Lyman breaks at 5.6 < z < 6.9, based on calculations 
from the population synt hesis model code GALAXEV 
( Bruzual fc Chariot] 120031) and intergalactic absorption 
( 




z'(mag) 

Fig. 1. — A color-magnitude plot of detected objects (dots) 
brighter than the 3<r limiting magnitude (z' = 27.09) in z'-band. 
The horizontal dotted line shows the color criterion i' — z' = 1.5 
for z ru 6 galaxy candidates, and the gray region shows i'-band 
magnitudes fainter than 1<t limiting magnitude (i' = 28.91). The 
open circles are the 258 z ~ 6 galaxy candidates. 

3 rd criterion was imposed only for objects detected in J- 
band. Only nine objects of 267 i'-dropout objects were 
detected in the J-band image, and eight objects with 
J > 0.8 were removed from the sample. The 4 th cri 



Madaulll995l) . In GALAXEV. we simulated a large va- terion was also required to keep the contamination rate 

low. One object, which did not satisfy the 4 th criterion, 
was removed from the sample. Finally, we obtained 258 
z ~ 6 galaxy candidates. 

We evaluated the contamination rate for these color- 
selection criteria. The sources of the majority of the 
contamination were M/L/T dwarfs and z ~ 1.3 old el- 
liptical galaxies, the latter being able to satisfy the 1 st 
criterion (i' — z 1 > 1.50) due to the 4000A Balmer break. 
We simulated old galaxy SEDs at z ~ 1.3 using the 
GALAXEV, assuming two relatively high metallicities 
(Z Q and 2.5Z Q ), and extracted model spectra with age of 
0.5—5 Gyr, applying reddening with E(B — V) of 0.5, 1.0, 
and 1.5. Based on these simulations, high amount of dust 
(E(B—V) > 1.5) is needed to meet the 1 st color criterion. 
According to the Extrem ely Red Objects (EROs) catalog 
bv lMivazaki et all (|2003t ). 96% of EROs are %' - z' < 1.5 
and no EROs have E(B- V ) > 1.5 at z - 1.3. Addition- 
ally, [MaIhotr^^F^L| ( 20051 ) show that i'-dropout objects 
with i' — z' > 1.3 do not include any EROs based on 
their spectroscopy. Furthermore, the 2 nd and 3 rd criteria 
arc supplementally used in order to discriminate low-z 
galaxies from z ~ 6 galaxies. Hence, we ignore the pos- 
sible contamination by low-z galaxies in the subsequent 
analyses and discussions. Secondly, we consider contam- 
ination by M/L/T dwarfs. According to lHawlev et al.l 



riety of galaxy spectral energy distributions (SEDs) using 
the Pad ova 1994 s imple stellar population model. We as- 
sumed a lSalpeterl ([1995) initial mass function with lower 
and upper mass cutoffs m^ =0.1 M Q and mu = 100 M Q , 
two metallicities (0.2 and 0.02 Z Q ), and instantaneous 
burst star formation history. We extracted model spec- 
tra with ages between 10 M yr and 740 M yr and applied 
the reddening law of lCalzetti et al.l ()2000D with E{B-V) 
of 0.00, 0.01 and 0.10. The 2 nd criterion was required be- 
cause the flux at the wavelength blueward of Lyu m it for 
z ~ 6 galaxies should be absorbed almost completely 
by the intergalactic medium. We did not require null 
detections in either the R- or i'-band because the con- 
tinuum blueward of Lya may not be totally extincted at 
z ~ 6 and because the Ly/3 and Ly7 emissions, if they are 
present in these galaxies, can fall in these bands. How- 
ever, it is interesting to note that the selected objects 
have fainter R magnitudes than the i' magnitude, or are 
not detected in the i?-band. The 3 rd criterion was re- 
quired because the spectrum of z ~ 6 LBG redward of 
Lya should be blue: UV-cont inuum slopes, j3 (f\ o c A* 3 ), 
are typically /3 < -2 (e.g.. iBouwens et al.1 12009). As 
mentioned in §2.1, our J-band image was > 2 mag shal- 
lower than the z'-band image; thus, almost no z ~ 6 
galaxies should be detected in the J-band image. The 
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Fig. 2. — The sky distribution of 258 z ~ 6 galaxy candidates, 
with surface number density contours. The z ~ 6 galaxy candi- 
dates are represented by filled circles whose size is proportional to 
the z'-band magnitudes (huge: 25.0 < z' < 25.5, large: 25.5 < 
z' < 26.0, medium: 26.0 < z' < 26.5, small: 26.5 < z' < 27.1). 
The lines correspond to contours of surface overdense significance 
from 6cr to Oct with a step of 2cr. North is up and east is to the 
left. The comoving scale projected to z = 6 is also shown along 
the axes. Gray regions are masked regions. The overdense region 
can be clearly seen at the southern edge of the plot. 



(2002), M/L/T dwarfs that meet the 1 st criterion should 
have a very red color [z 1 — J > 2). Thus, M/L/T dwarfs 
can be, in principle, detected in the J-band image such 
that they would be removed from our z ~ 6 galaxy can- 
didate sample with the 3 rd criterion. However, since the 
J-band limiting magnitude is > 2 mag shallower than 
that of the z'-band, in practice the 3 rd criterion is effec- 
tive to remove only the brightest or very reddest M/L/T 
dwarfs from our high-redshift sample. The 4 th criterion 
was also imposed to restrain the contamination rate by 
M/L/T dwarfs because the contamination rate is high at 
brighter than z' ~ 25. From th e star count model devel- 
oped by iNakaiima et aLl (|2000[ ). the contamination rate 
at z' > 25 is only about 6% at the galactic coordinates 
corresponding to our survey area. Therefore, we assumed 
a contamination rate in our z ~ 6 galaxy selection of up 
to a few percent, mainly consisting of contamination from 
fainter or bluer M/L/T dwarfs. 

In addition to these intrinsic contaminates, photomet- 
ric noise may scatter lower- redshift objects to satisfy 
our sample selection criteria. We performed the fol- 
lowing simple simula tion, the same as is perfor med in 
IWilkins et aLl (|2011al ) and lBouwens et all (|2011al ). to es- 
timate the contamination rate due to the photometric 
noise. We first randomly choose brighter objects at 
22.0 < z' < 22.5, dim these bright objects so as to match 
the magnitude distribution to those of our i'-dropout at 
25.0 < z' < 27.09 objects by scaling the flux, then dis- 



TABLE 2 

Photometric Properties of the 15 Spectroscopic 

Confirmed Galaxies 



ID 



RA 
(J2000) 



DEC 
(J2000) 



(mag) 



(mag) 



l a 


13:24:31.8 


+27:18:44.2 


25.90 + 0.04 


1.48 


2 


13:24:18.4 


+27:16:32.6 


25.69 + 0.03 


1.70 


3 


13:24:25.2 


+27:16:12.2 


27.02 + 0.09 


> 1.89 


1 


13:24:30.2 


+27:14:13.5 


26.81+0.08 


1.96 


5 b 


13:24:26.0 


+27:16:03.0 


26.50 + 0.06 


1.71 


6 


13:24:21.3 


+27:13:04.8 


25.91+0.04 


2.65 


7 


13:24:29.0 


+27:19:18.0 


26.50 + 0.06 


1.60 


ga 


13:24:28.1 


+27:19:32.7 


26.10 + 0.04 


1.54 


9 


13:24:26.5 


+27:15:59.7 


25.47 ± 0.03 


1.91 


10 a 


13:24:31.5 


+27:15:08.8 


25.91 + 0.06 


1.95 


11 


13:24:26.1 


+27:18:40.5 


26.59 + 0.06 


> 2.32 


12 


13:24:31.6 


+27:19:58.2 


26.47 ± 0.06 


2.07 


13 


13:24:44.3 


+27:19:50.0 


26.43 ± 0.06 


> 2.49 


14 c 


13:24:20.6 


+27:16:40.5 


27.01 + 0.09 


> 1.90 


15 a 


13:24:32.6 


+27:19:04.0 


25.54 + 0.03 


1.47 



a These do not satisfy all of our selection criteria. 

The same object as No. 6 in Table 1 of Jiang; et al. (2011). 
c The same object as No. 7 in Table 1 of Jiang et al. (2011). 



tribute these artificially dimmed objects on the original 
image as we did in §2.1. In this simulation, we simply as- 
sumed that the intrinsic color distribution of faint objects 
is almost the same as that of brighter objects whose pho- 
tometric noise should be negligible. The bright sources 
are chosen from only slightly brighter magnitude ranges, 
in which they have higher S/N but are expected to have 
similar color distribution with our i'-dropout sample. We 
extracted the artificial objects using SExtractor and im- 
posed our color criteria of the i'-dropout objects. The 
number of artificial objects in a given z'-magnitude in- 
terval was chosen to be the same as the observed number 
of object in the same z'-magnitude interval. We finally 
found that 39 artificial objects that met our color criteria, 
i.e., 15% of 258 i'-dropout objects may be contaminated 
due to the photometric noise. It should be noted that 
these contaminants appeared almost randomly over the 
image in the repeated simulation, although the depth of 
J-band image is not uniform over the field. This suggests 
that the contamination rate is almost homogeneous over 
the survey field and did not change the overdensity sig- 
nificance estimated in the next section. 

3. SKY DISTRIBUTION AND OVERDENSITY 

SIGNIFICANCE 

Figure [5] shows the sky distribution of the 258 z ~ 6 
galaxy candidates. We found an apparently overdense 
region, centered in the southern part of the field. To 
determine the overdensity significance quantitatively, we 
estimated the local surface number density by counting 
z ~ 6 galaxy candidates within a fixed aperture of 2.'1 
radius, corresponding to 5Mpc at z = 6. These aper- 
tures were distributed over the SDF in a grid pattern at 
intervals of 20 arcsec. We assumed that the local sur- 
face number density in masked regions was the same as 
the mean surface number density. The mean and the 
dispersion, a, of the number of z ~ 6 galaxy candi- 
dates in a circle was found to be 4.0 and 2.6 per cell, 
respectively. Based of these values, we also plotted sur- 
face number density contours in Figure [2j We found 
that our region of interest appeared overdense at the 6<7 
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Fig. 3. — Spectra of eight protocluster member galaxies. The 
vertical dashed lines indicate the peak of Lya emission line. Blue 
solid lines represent sky lines. The numbers at the upper left corner 
are the IDs (Column 1 of Table 3). 



significance level beyond the mean surface number den- 
sity at the peak. We designated an "overdense region" 
6' x 6' (14 x 14 Mpc 2 ) centered on the peak of overdensity. 
This encompasses a region with an overdensity of more 
than 3cr significance. This overdense region included 30 
z ~ 6 galaxy candidates. The size of the overdense re- 
gion was al most the same as that of the protoc l usters 
at z = 2-5 (jVenemans et al.l 12007( 1 . lOuchi et al.l ((20051 ) 
found two protoclusters whose sizes were ~ 7 Mpc di- 
ameter at z = 5.7 i n the Subaru/XMM-Newton Deep 
Survey. iKang fc iinl ((20091) also found an overdense re- 
gion with ~ 10 Mpc diameter at z = 3.7 in Chandra 
Deep Field-South. From these comparisons, the size of 
our overdense region is almost consistent with previous 
observational results of protoclusters at lower-z. 

We also co mpared our overdense region with theoret- 
ical models. lOverzier et al.l ((2009a( l constructed a very 
large (~ 4° x 4°) mock redshift survey of star- forming 
galaxies at z ~ 6, using the Millenn ium Run dark mat- 
ter simulation (jSpringel et al.l (2005) and semi-analytic 
prescriptions of galaxy formation. In this mock survey, 
the limiting magnitude of z'-band was 26.5 mag. They 
predicted a typical protocluster diameter of not more 
than ~ 25 Mpc at z ~ 6. The size of our overdense 
region (< 14 Mpc) is consistent with this theoretical re- 
sult. Furthermore, they calculated the probabilities of 
finding overdense regions and the fraction of i'-dropout 
objects marked as protocluster member galaxies. Ap- 
plying the r'-dropout selection criteria of Overzier et al. 
(2009a) to our catalog, the probability of finding a region 
with an overdensity significance higher than ours within 
a 3.4 x 3.4 arcmin 2 (the FOV of HST/ACS) is only about 
0.3%; the local sur face number density wa s estimated fol- 
lowing the way of lOverzier et al.l ((200 9a). Furthermore, 
it is estimated from these results that 50% of z ~ 6 pho- 
tometric candidates in this overdense region belong to the 
protocluster, though the scatter is considerable (lc range 
is between 0.3 and 0.7). The overdense region contains 
fainter galaxies at 26.5 < z' < 27.1; the probability of 
these being real protocluster members is not predicted. 
The discovery of an overdense region with such a high 
statistical significance would be extremely rare, and half 
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Fig. 4. — Same as Figure[3] but for the seven non-member galax- 
ies. 

of all galaxies in the region are expected to be protoclus- 
ter galaxies. Therefore, there is a high probability that 
the overdense region includes a protocluster at z ~ 6. 

4. SPECTROSCOPIC CONFIRMATION 

Our discovery of an overdense region might be at- 
tributed to a mere chance of alignment, because the 
dropout technique samples a broader range of redshifts. 
It might also be an incidental result of highly clustered 
contaminating populations. To investigate this, we car- 
ried out follow- up spectroscopy on the overdense region, 
using FOCAS (Kashik awa et all 120021) on the Subaru 
telescope in Multi-Object Spectroscopy (MOS) mode. 
The details of this observation will appear elsewhere 
((Shibuvaet al.l 12011( 1. The FOCAS observations were 
made with a 555 mm -1 grating and an 058 order-cut fil- 
ter, giving spectral coverage of 7500-10450A with a res- 
olution of 0.74 A pix -1 . The 0'.' 8- wide slit gave a spec- 
troscopic resolution of 5.7A (R~1500). We targeted 27 
z ~ 6 galaxy candidates in the overdense region. Addi- 
tionally, we allocated slits to eight objects that did not 
meet our photometric criteria perfectly, i.e., i'—z' > 1.40, 
z' < 25.0, or B-band magnitude slightly brighter than 
the 2cr limiting magnitude, to evaluate our selection cri- 
teria. The data were obtained on 2010 March 19 and 
20. Sky conditions were good with a seeing size of 0'.'5- 
l'.'O. We obtained eight exposures, for a total integration 
time of four hours. Between exposures, the telescope 
was dithered along the slit to enable more accurate sky 
subtraction. Long slit exposures of standard star HZ 44 
were used for flux calibration. The data were reduced in 
a standard manner with IRAF. 

We carefully discriminated real emission lines from sky 
lines by examining 2D and ID spectra, and detected sin- 
gle emission lines from 10 objects of 35 targets. Our spec- 
troscopic observations were carried out with R ~ 1500, 
which is a high enough resolving power, in principle, 
to distinguish single Lya emission from [On] doublet 
(AA = 6.3 A at ~ 8500 A), although it is difficult in prac- 
tice to discriminate between the features if the emission 
is very faint. Rather than relying on the resolving power 
of our spectroscopy, we use the emission line profile as 
the primary diagnostic for distinguishing Lya at high 
redshift from [On] emission. We calculated the weighted 
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Fig. 5. — Thumbnail images (~ 15" X 15") of eight protocluster 
member galaxies in six passbands. From left to right, the B-, V-, 
R-, i'-, z'-, and J-band images are shown. North is up and east 
is to the left. The numbers at the upper left corner are the IDs 
(Column 1 of Table 3). 



skewness, S w (Kashika wa et all 120061 ) . which is a good 
indicator of the Lya emission-line asymmetry. The re- 
sults are listed in Table 3. We found that most of the 
detected emission lines had S w > 3, indicating that they 
are certainly Lya emission lines at z ~ 6; however, some 
emission lines did not meet the Lya criterion. These 
line profiles may be affected by strong sky lines in this 
wavelength range, or low S/N ratio data could prevent 
accurate determination of the skewness. As mentioned 
in §2.1, it is unlikely that [On] emitters at z ~ 1.3 sat- 
isfy the photometric selection criterion of i' — z' > 1.50. 
Although we cannot completely deny the possibility that 
objects with S w < 3 were [On] emission lines, the single 
emission lines from our i'-dropout sample were plausibly 
identified as Lya emission at z ~ 6. We hereafter regard 
the single emission line detected from our i'-dropout sam- 
ple as Lya emission. 

We detected Lya emissions from seven objects of 27 
targeted z ~ 6 galaxy candidates that met all of our 
photometric criteria. We did not detect any signal from 
the remaining 20 objects. No apparent contaminations, 
such as M/L/T dwarfs or low-z galaxies, were detected. 
We allocated slits to eight candidates that did not meet 
our photometric criteria: three of these were identified 
as Lya emitters at z ~ 6, three were apparent stars. We 
did not detect any signal from the remaining two objects. 

We also obtained five Lya emission lines at z ~ 6 in the 
overdense region from previous observations , carried out 
with Subaru/FOCAS or Keck/DEIMOS (jFaber et all 
|2003l ). These spectra were obtained through pre- 
vious MOS observations on several different projects 
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Fig. 6. — Same as Figure(5] but for the seven non-member galax- 
ies. 

(iNagao et all 120071: lOta et all 120081: Kashikawa et al.1 



l201lll,Tiang et all 1201 lHShibuvaet al1l2011l ). Our z ~ 6 
galaxy photometric candidate sample includes four of 
these five galaxies. The remaining z ~ 6 galaxy was 
so close to a bright star that we could not detect it. 

In sum, we spectroscopically confirmed 15 z ~ 6 galax- 
ies, including four serendipitous objects, in the overdense 
region. The details of all 15 z ~ 6 galaxies are provided 
in Tableland 3. The spectra and thumbnail images of 
all 15 z ~ 6 galaxies are shown in Figure |31 HI and [SJ 
[HI respectively. All the z ~ 6 galaxies with measured 
redshifts are identified by their Lya emissions. 

The spectroscopically undetected targets were ex- 
pected to have a Lya flux that is too low to be detected 
(or non-existent). Our sample was selected based on the 
UV continuum flux using broad-band images. Therefore, 
these were not necessarily required to have Lya emission. 
As noted in §2.1, our z'-band image, with a 30-hour-long 
integration time, was stacked from all the data taken for 
8 years. Therefore, there is little chance that they are 
variable objects, such as supernovae or active galactic 
nuclei that incidentally increased their luminosity dur- 
ing the z'-band observation . The misalignment of slits 
is also quite unlikely because we used the MOS mode 
with enough alignment stars, and we actually succeed in 
detecting signals in some of slits in the same MOS mask. 
The detection rate of Lya emission of i'-dropout candi- 
dates is 75% at z' < 26, 38% at 26 < z' < 26.5, and 22% 
at z' > 26.5. 

We made a composite spectrum of all 15 Lya emission 
lines. Each galaxy spectrum was shifted to its rest-frame 
and normalized using the flux at peak. Then, we created 
a weighted average of all spectra in which the weight was 
determined from the S/N ratio. The resulting compos- 
ite spectrum is shown in Figure clearly revealing an 
asymmetric profile, a characteristic feature of Lya emis- 
sion lines at high redshift. The weighted skewness was 
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Fig. 7. — The composite spectra of spectroscopically confirmed 
galaxies. The blue, red, and green lines indicate that of all 15 galax- 
ies, eight protocluster member galaxies, and seven protocluster 
non-member galaxies, respectively. The red-shaded region shows 
the la variance on the composite spectrum of member galaxies. 

found to be S w = 8.76 ± 1.05. 

The redshifts of eight galaxies (ID=3-10 in Table 3) of 
15 LAEs in the overdense region lie between z = 5.984 
and z = 6.027, apparently concentrating at z = 6.01 
(Az < 0.05). Figure [8] shows the redshift distribution of 
15 galaxies. The redshift concentration at z — 6.01 is 
about seven times higher than the number density ex- 
pected from a uniform distribution. Moreover, it should 
be noted that all eight galaxies clustered in redshift have 
S w > 3, suggesting that these are plausible Lya emitters 
in a protocluster at z ~ 6. Figure |9] shows the three- 
dimensional distribution of these eight galaxies. These 
eight galaxies are strongly clustered, with separations 
less than 16 Mpc. Thus, this region is almost certain to 
be a galaxy protocluster at z = 6.01. If true, this proto- 
cluster would be the highest redshift large scale structure 
probed by the combination of the i'-dropout technique 
and the unique wide-field imaging capability of Suprime- 
Cam. 

5. DISCUSSION 

5.1. Uncertainties in the Overdense Significance 

We applied a color threshold i' — z' > 1.50 to select 
i'-dropout objects, although the critical threshold does 
not perfectly eliminate other contaminations, and photo- 
metric errors dilute the discrimination. We evaluated the 
effect of the critical color threshold on the significance of 
the overdense region by changing the threshold of the 1 st 
criterion from 1.5 to 1.3, 1.4, 1.6 and 1.7 in the sample se- 
lection. We re-estimated the significance of the overdense 
region relative to spatially uniform distribution with the 
same procedure as in §3. The overdense region was al- 
ways uniquely identified as the highest-density region in 
the SDF, and the peak significance fluctuated between 
5a and 7a, indicating that the significance was statisti- 
cally robust. As noted in §2.2, the contamination from 
M/L/T dwarfs could be varied from field to field due to 
non-uniform J-band depth over the field. However, since 
the overdensity region is mostly located in the deepest 
portion of the mosaic J-band image, the overdensity sig- 
nificance should not be incidentally enhanced by the con- 
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Fig. 8. — The redshift distribution of 15 spectroscopically con- 
firmed galaxies. The bin size is Az = 0.05. The solid (blue) line 
shows the selection function of our i'-dropout selection assuming a 
uniform distribution normalized to the total number of confirmed 
emitters. The inset is a close-up of the protocluster redshift range, 
with a bin size of 400kms — 1 . The horizontal axis, Av, represents 
radial velocity. The eight protocluster members are distributed be- 
tween the vertical dashed (blue) lines. The dotted lines show the 
location of sky lines if redshifts are determined by Ly« emission. 



tamination. 

The spectroscopic identification of Lya emitters has 
some uncertainty. In §4, we regarded all single emis- 
sion lines detected from our i'-dropout sample as Lya 
emission, although we found some emissions with asym- 
metries which were not high enough to be conclusively 
identified as Lya. The objects labeled as ID=1 and 15 in 
Table 3, have S w < 3. These might be [On] emitters. It 
should be noted that all eight protocluster galaxies have 
S w > 3. Even if the objects with S w < 3 were deleted 
from the z ~ 6 galaxy, the concentration of the eight 
galaxies still remains. Furthermore, the detection limit 
of Lya emission is not, generally, a uniform function of 
wavelength. This is due to the effect of strong night sky 
lines in certain spectral window, which can overwhelm 
the faint Lya emission we are attempting to detect. The 
fact that all of our protocluster members are detected 
within 8520 ± 30 A, a window not free of bright airglow 
lines, further bolsters our confidence that the redshift 
clustering of our targets has an astrophysical origin and 
is not a consequence of selection effects resulting from 
night sky emission. 

5.2. Protocluster Properties 

We estimated the observed properties of the spectro- 
scopically confirmed galaxies, properties such as the Lya 
fluxes, the UV continuum magnitudes, and the Lya rest- 
frame equivalent widths (EWs). Although faint contin- 
uum flux could not be detected in the observed spectra, 
the UV continuum magnitudes were estimated from the 
z'-band magnitudes by subtracting the spectroscopically 
measured Lya fluxes and assuming flat UV continuum 
spectra (/„ = constant). The 15 z ~ 6 galaxies are di- 
vided into two groups; the protocluster "members," in- 
cluding ID=3-10, and "non-members," including ID=1, 
2, and 11-15. We compared the UV continuum luminos- 
ity, Lya luminosity, rest-frame EW and FWHM of these 
two groups. The average of these parameters between the 
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Fig. 9. — Three-dimensional distribution of protocluster member galaxies. The upper left panel shows the spectroscopically observed 
region. The lines and points are the same as those in Figure [2] and the crosses and triangles are spectroscopically observed objects, the 
Lya emissions of which were detected and not detected, respectively. The numbers near the points are the IDs (Column 1 of Table 3). 
The lower right panel is a closeup of the central region of the protocluster. The vertical axis in the lower left panel and the horizontal axis 
in the upper right panel, A, represent radial distance with an origin at z = 6. Note that in this figure, we regarded redshift difference as 
distance difference, although redshift difference might be caused by peculiar velocities inside the protocluster. 



members and the non-members agree with each other 
to within la scatter. The composite spectrum of the 
members is shown in Figure [71 and its weighted skew- 
ness is S w — 8.86 ± 2.14, whereas the weighted skew- 
ness of the composite spectrum of the non-members is 
S w = 8.70 ± 2.04. Thus, the properties of the members 
are not found to be significantly different from those of 
non-members despite the difference in the environments 
between the two samples. However, we should note that 
the scatter in each of these properties is too large to con- 
clude definitely and that additional spectroscopic con- 
firmations are required for more accurate comparisons. 
It is also worth noting that we have not found appar- 
ent signatures of AGN activity for any of the member 
galaxies. At the center of the protocluster, there is a 
galaxy (ID=9) whose z'-band magnitude of z' — 25.47 is 
> 0.9 mag brighter than the average of all members. This 
bright galaxy has an SFR(UV) = 110.4 Meyr -1 , which 
seems too small for strong AGN activity. We also did not 
detect any N V A1240A emission from the object. This 
was the only accessible line that could have indicated 
AGN activity in our optical spectrum, although deeper 
spectroscopy is required to detect its AGN signature, if 



present. In §2.2, we mentioned one object that was too 
bright to meet the 4 th criterion. The brightness of this 
object indicates that it may be an AGN, but it is not 
located in the protocluster region. The excess of AGNs 
or strong AGN activity were found in low redshift pro- 
toclusters, possibly resulting from merging, which would 
be enhanced in such a high-density e nvironment (e.g., 
iLehmer et all 120091: iKuiper et all 1201 1[) . Unlike the low 
redshift protoclusters, strong AGN activity could not be 
found in our protocluster region, implying that the mem- 
ber galaxies of this z ~ 6 protocluster may be too young 
to have the requisite conditions needed for significant 
AGN activity. 

5.3. The Mass Estimate of the Protocluster 

The eight clustering galaxies are spread over a 3f5 x 6.'5 
region, corresponding to 8.3 x 15.5 Mpc 2 . The radial 
velocity dispersion, estimated simply from the disper- 
sion of measured wavelength of the Lya emission line, 
is ay = 647 ± 124kms~ 1 . The error was estimated 
from the uncertainty of the wavelength of the Lya emis- 
sion line, taking into account the wavelength resolu- 
tion and the strong sky lines. Furthermore, the res- 
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onant scattering of Ly a line could lend to additional 
error. iVenemans et al.l (|2007f) found that the velocity 
dispersions of protoclusters at z ~ 2—3 and 4—5 are 
o> ~ 500—1000 and 300 kms -1 , respectively. The dark 
matter velocity dispersion is predicted to increase with 
cosmic time, ba sed on a numeri cal simulation of cluster 
evolution (e.g., lEke et ail [T998) . This is roughly con- 
sistent with observations at z = 2—3 and 4—5, and is 
predicted to reach to a r ~ 200 km s^ 1 at z ~ 6. Our 
estimate of the velocity dispersion at z ~ 6 is much 
larger than this. We will discuss this discrepancy in 
the next subsection. The radius of the protocluster is 
r ~ 1 Mpc in physical units. Thus, the virial mass, 
Mvi r ~ 3cr^r/G, where G is the gravitational constant 
and r is the radius of the structure in physical units, 
would be 2.9JIj'o x 10 14 M Q , simply assuming that the 
protocluster region was in virial equilibrium. However, 
the protocluster is expected to be far from virial equilib- 
rium. 

Instead, we used another method to estimate the 
mass , specifically, M =_pV(l + S m ) ([Steidel et al.l 
119981: IVenemans et al.l 120051) . where p is the current 
mean matter density of the universe, V is the comov- 
ing volume of the protocluster and S m is the mass 
overdensity of the protocluster. The mean density is 
p = 4.1 x 10 10 Mq Mpc -3 for the cosmological param- 
eters assumed in this paper. The mass overdensity, 
S m , is related to the galaxy overdensity 5 sa x, through 
1 + b8 m — C(l + <5 g ai), where b is the bias parame- 
ter and C takes into accou nt the redshift space dis- 
tortions dSteidel et al.l I1998D. Semi-ana lytical models 
(jBaugh et al.lll9 99; Kauffmann et al. 199 9*) and hydrody- 
nami c simulations (Blanton et al. 2000; Yoshikawa et al. 
120011 ) predict that the bias parameter will increase with 
increasing redshift. The bias paramet er is predi c ted to 
be b = 4.0 ± 0.8 for z = 5.7 LAEs bv lOrsi et all (f2008h 
based on their semi-analytical model in the ACDM uni- 
verse. Thus, we assumed the bias parameter for LAEs 
at z ~ 6 to be b — 4.0 ± 0.8. It should be noted that 
the observe d bias parameter o f LAEs still has a large 
uncertainty. IQuchi et al.l ()2003D derived the bias param- 
eter t o be 6 = 5.0 for z ~ 4.8 LAEs, and iKovac et al.l 
(|2007f) derived b = 3.7 for z - 4.5 LAEs, so these re- 
sults should be interpreted with caution. C can be ap- 
proximated by C = 1 + / - /(l + <5 m ) 1/3 (ISteidel et al.l 
1998), where / is the rate of growth of perturbations 
at the redshift of the protocluster. At high redshift, 
this value is close to 1 (La hav et al.lll99l"f) . If we adopt 
the volume of the protocluster as V — 2.3 x 10 3 Mpc 3 , 
with six protocluster galaxies being identified in the vol- 
ume (this excludes ID=8, 10, which are the serendipi- 
tous identifications), the protocluster number density is 
^cluster = (2.6 ± 1.1) x 10~ 3 Mpc -3 . Using the mean sur- 
face number density (§3.1), the field number density is 
nfidd = (1.5 ± 0.1) x 10~ 4 Mpc -3 . These yield a galaxy 
overdensity of <5 gal = n c i ustor /n fio id - 1 = 16 ± 7, corre- 
sponding to 6 m = 2.1 ± 0.9. Therefore, the protocluster 
mass is M = (2.9 ± 0.8) x 10 14 M Q . 

Although it is interesting that the masses estimated 
from two different methods are consistent with each 
other, each mass estimate has large uncertainties. In 
the first method, we assumed virial equilibrium, which 
is almost certainly not true for such a structure in the 



early phase of its evolution. The second method re- 
tains a large uncertainty in estimating mass overdensity 
from the galaxy number density. Additionally, we re- 
garded differences in redshift as difference in distance 
within the protocluster. However, clustering structure 
may break the Hubble flow; in this case, the proto- 
cluster volume would be overestimated. Nevertheless, 
based on these mass estimations, the halo mass of the 
present-day descendant can be predicted in a statisti- 
cal sense using the extended Prcss-Schechter (EPS) for- 
malism (jBond et all 119911: IBowerlll99"i"l) . We find that 
the mass of ~ 2.9 x 10 14 M Q at z = 6 will become the 
mass of (3.5-5.8) x 10 15 M Q (68% range of distribution 
function) at z = 0. The predicted mass of the descen- 
dant is higher than the typical mass of local rich clus- 
ters ~ 9 x 10 14 M Q and comparabl e to t he most massive 
cluster - 4 x 1O 15 M (|Wen et al.l 120101 ): thus, the high 
concentrations at z ~ 6 might become one of the richest 
clusters in the local u niverse. The theoretical work by 
IQverzier et al.l (|2009al) predicted that the most massive 
halo masses are ~ 10 13 M Q at z ~ 6, and these evolved to 
rich clusters with halo masses of ~ 10 14 — 10 15 M in the 
local universe. The derived mass, which is overestimated 
compared with the prediction, would imply two ideas: 
either the protocluster is in a very early stage of cluster 
formation, which leads such a large velocity dispersion, 
or it is a progenitor of an extremely massive cluster in 
the local universe. 

5.4. Protocluster Structure 

The estimated velocity dispersion of the protocluster 
is much larger than those at lower redshifts, and larger 
than the dispersion predicted by CDM models. This 
discrepancy could be due to the distinguishing three- 
dimensional structure of the protocluster. Figure[9]shows 
the three-dimensional galaxy distribution of the proto- 
cluster, with possible substructures. The spectroscopi- 
cally confirmed galaxies tend to be distributed on the 
surface of a sphere ~ 8 Mpc in radius, but a void con- 
taining few spectroscopically identified galaxies seems to 
be present at the center of the sphere. The z ~ 6 galaxy 
candidates in the void are generally faint in z'-band of 
z' > 26.5, in which the detection rate of Lya emission 
is only 22%, as described in §4. In contrast, the frac- 
tion of Lya emitters among LBGs at the faint end of the 
LBG luminosity function is gener ally higher than that of 
the bright end (jStark et alll201lD . Although this result 
seems to contradict ours, Stark et al. made much deeper 
spectroscopic observations than ours. In our bright sam- 
ple with z'-band magnitude (25.0 < z' < 26.5), the LAE 
fraction with EW > 25A is 30±17%. Th i s is co mparable, 
within the error, to that of Sta rk et all (1201 ip. which is 
20±8.1%. If the same trend of the fraction in lStark et all 
(|2011l ) applied to our faint sample, about half of the spec- 
troscopically undetected galaxies with z' > 26.5 are ex- 
pected to have Lya emissions with small EW (< 85 A), 
which is difficult to detect in our observations. There- 
fore, it is likely that this effect is the source of the void 
seen at the center of the protocluster. 

Although it is not clear whether the z ~ 6 galaxy can- 
didates in the void region, including spectroscopically 
undetected galaxies, are protocluster members or non- 
members, we now consider two extreme scenarios. First, 
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we suppose that all z ~ 6 galaxy candidates in the void 
region are real protocluster members, forming a large 
protocluster with a radius of ~ 8 Mpc. Second, we sup- 
pose that none of them are protocluster members, that 
the eight protocluster galaxies are divided into small sub- 
groups. In the void region, we found 17 z ~ 6 galaxy can- 
didates without Lya emission lines. If their redshift dis- 
tribution is similar to that of our spectroscopically con- 
firmed samples (the histogram of Figure|8]), eight galaxies 
would be expected to be protocluster members. The to- 
tal number of protocluster members is 16 in this case, 
and this supports the first picture. In contrast, if the 
redshift distribution of the void galaxies mimics a uni- 
form distribution (the line of Figure [8]) , only one galaxy 
is expected to be a protocluster member galaxy. In this 
case, we would assume the second picture was accurate. 
As discussed in §4, the fraction of the real protocluster 
members in the overdense region is expected to be 0.5. 
We cannot determine which hypothesis is likely bas ed on 
the numerical simulation by lOverzier et al.l (|2009a|) . be- 
cause the fraction of protocluster member galaxies has a 
large scatter. 

In the first case, the total number of protocluster mem- 
ber galaxies would be 25 at most. Galaxies located in the 
void are considered to be LBGs, either with weak Lya 
emission, or with no Lya emission. In contrast, LBGs 
with strong Lya emission are predominantly distributed 
in the outer region of the protocluster. In this case, our 
estimate of velocity dispersion of 647kms~ 1 , using only 
outer protocluster galaxies, is an overestimate; the ac- 
tual velocity dispersion would likely be m uch smaller. 
LBG s are generally older than LAEs ( cf., lEvles et al.l 
[20071 IStarket al.l[200l IQno et al.ll2010h . If this is the 
case in this protocluster, inner galaxies in the protoclus- 
ter may have evolved faster than the outer galaxies, and 
may have already finished their early LAE phase. Based 
on this picture, the radial dependence of galaxies' prop- 
erties has already appeared even in the z ~ 6 proto- 
cluster. While these results are suggestive, it would be 
interesting to obtain deeper spectroscopic observation of 
such systems to definitively see whether such an envi- 
ronmental effect is already in place at such early epochs. 
A similar feature was found in a protocluster at z ~ 5 
(Kashikaw a et al.ll2007fl , though it was not distinguished 
from the possible negative feedback effect from the QSO 
at the center of the protocluster. 

In the second case, spectroscopically confirmed galax- 
ies distributed on the surface of a sphere seem to split 
into three subgroups: the 1 st , containing ID=3-6; the 
2 nd , containing ID=7 and 8; and the 3 rd , containing 
ID=9 and 10. These subgroups might be in the pro- 
cess of forming a larger protocluster through merging. 
In this case, the overall protocluster system is expected 
to be far from virial equilibrium; thus, the dynamical 
mass obtained from the velocity dispersion of eight LAEs 
was a significant overestimate. The velocity dispersion 
of a subgroup including only four galaxies (ID=3-6 in 
Table 3) is estimated to be 224 ± 152 km s _1 , which 
is in a good agreement with the CDM model predic- 
tion at z ~ 6. The virial mass of the subgroup is 
2.0l^'g x 10 13 M Q , and the mass of the descendant is 
predicte d to be (3.0-8.9) x 10 14 M Q by the EPS for- 
malism. lOuchi et al.l (|2005) found two protoclusters at 



z ~ 5.7, and their velocity dispersions were ~ 180 and 
150kms _1 . We may be looking at a site of cluster for- 
mation prior to collapse, in which three small groups of 
galaxies begin to merge with each other. Such a merger 
may trigger AGN activity at the cente r of the protoclus- 
ter in the future. iKuiper et al.l (|2011h found a merging 
protocluster at z ~ 2.2, and suggested that a merger 
could enhance star formation and AGN activity. Either 
scenario presents intriguing possibility. 

We regarded the redshift difference as the difference 
in radial distance when discussing two scenarios above; 
however, these differences might be caused by the pe- 
culiar velocity inside the protocluster. More detailed 
investigations are necessary to reveal primitive cluster 
formation. A series of redshift determinations for the 
galaxies in the void, derived from deep spectroscopy, is 
the most straightforward approach for distinguishing be- 
tween these two ideas. If the first case is true, galaxies 
in the void region would have older ages, higher stellar 
masses, higher metallicities, and larger amounts of dust 
extinction. We could confirm those points with observa- 
tions taken at longer wavelengths to trace the rest-optical 
frame of z ~ 6 galaxies, followed by accurate SED fitting. 
We plan to address this issue in detail in the future. If 
the second case is true, it would be interesting to in- 
vestigate the possibility of additional subgroups further 
outside the protocluster. Spectroscopic observations de- 
signed to locate outer galaxies should reveal additional 
subgroups or filamentary structures, closely related to 
the initial large-scale structure of the universe. 

6. CONCLUSIONS 

In this paper, we have presented the discovery of a 
protocluster at z ~ 6 in the SDF based on eight spec- 
troscopically confirmed member galaxies. This discovery 
was achieved in a "random field" using wide-field imag- 
ing; that is, QSOs or RGs were not used as protocluster 
probes. The z'-dropout surface number density of the 
protocluster region is overdense at the 6<r level compared 
with the mean density of the whole SDF. The size of the 
overdense region is ~ 6' x 6' (14 Mpc x 14 Mpc). The 
protocluster galaxies are distributed in a narrow redshift 
range of Az < 0.05 around z = 6.01. The concentra- 
tion in redshift space is seven times higher than that 
expected from a uniform distribution. This region is cer- 
tainly a protocluster at z = 6.01, the highest redshift 
such a structure has ever been observed. The differences 
in galaxy properties, such as Lya luminosity and UV con- 
tinuum magnitude, between protocluster members and 
non-members (i.e., field galaxies) at z ~ 6 cannot be de- 
tected, although the small sample size still prevents us 
from forming a definite conclusion. 

The estimated velocity dispersion of the protocluster is 
a r = 647 ± 124 kms" 1 . Assuming virial equilibrium, the 
protocluster mass is estimated to be 2.9JZ^q x 1O 14 M . 
Based on the EPS formalism, we expect that, if viri- 
alized, this protocluster will become one of the richest 
clusters (~ 4.7 x 10 15 M Q ) in the local universe. From 
the three-dimensional distribution of the eight protoclus- 
ter galaxies, we proposed two different scenarios for the 
protocluster structure: either the protocluster is already 
mature, with old galaxies at the center, or it is still imma- 
ture and composed of three subgroups merging to become 
a larger cluster. Further deep spectroscopic observations 
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or multi-wavelength imaging to trace the rest-frame op- 
tical wavelengths would provide a clearer picture of the 
protocluster structure, as well as providing constraints 
for the galaxy /stellar populations of protocluster mem- 
bers. A statistical sample of z ~ 6 protoclusters based on 
wider surveys would allow a direct comparison of the ob- 
served protocluster abundance in the early universe with 
the predictions of the ACDM model. 

This paper is based on data collected with the Subaru 



Telescope, which is operated by the National Astronomi- 
cal Observatory of Japan. We are grateful to the Subaru 
Observatory staff for their help with the observations. 
We thank the anonymous referee for valuable comments 
and suggestions which improved the manuscript. This re- 
search was supported by Japan Society for the Promotion 
of Science through Grant-in- Aid for Scientific Research 
23340050. 

Facilities: Subaru (SuprimeCam, FOCAS), Keck:II 
(DEIMOS), UKIRT (WFCAM) 
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